Ground states of unfrustrated spin Hamiltonians satisfy an area law 
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We show that ground states of unfrustrated quantum spin- 1/2 systems on general lattices satisfy an entan- 
glement area law, provided that the Hamiltonian can be decomposed into nearest-neighbor interaction terms 
which have entangled excited states. The ground state manifold can be efficiently described as the image of a 
low-dimensional subspace of low Schmidt measure, under an efficiently contractible tree-tensor network. This 
structure gives rise to the possibility of efficiently simulating the complete ground space (which is in general 
degenerate). We briefly discuss "non-generic" cases, including highly degenerate interactions with product 
eigenbases, using a relationship to percolation theory. We finally assess the possibility of using such tree tensor 
networks to simulate almost frustration-free spin models. 



I. INTRODUCTION 

An important insight in the study of quantum many-body 
systems is related to the observation that common states that 
naturally occur do not quite exhaust the entire Hilbert space 
available to them, but instead a much smaller subspace. This 
insight is at the heart of powerful numerical methods that have 
been devised in recent years. Ideas such as the density-matrix 
renormalization group approach, and new ideas that allow for 
the simulation of higher-dimensional quantum lattice models 
[OJO], work exactly because they model well quantum states 
that in a certain sense have little entanglement. More pre- 
cisely, the states which are tractible by these approaches sat- 
isfy what is called an area law ifTl l4l fI71 . so the entropy of 
a subregion scales at most as the boundary area of that re- 
gion (for a review, see Ref. (T)). For practical purposes, and 
in particular for ID systems, these methods in particular give 
accurate accounts of ground state properties. 

Now, not all ground states of local quantum lattice mod- 
els can be efficiently approximated. This holds even true for 
ID chains: indeed, one can construct models for which ap- 
proximating the ground state energy is provably NP-hard lfT8l 
— albeit using a fairly sophisticated construction involinvg 
large local dimensions 1201 . An important feature of these 
constructions is that the difficulty of their solution appears to 
be strongly related to whether the system is frustrated or not. 
This suggests that whether or not the system is frustrated is an- 
other criterion for whether a quantum lattice model should be 
considered "easy" or "hard", in addition to its ground states 
having "a lot" or "little" entanglement. This intuition that 
frustrated systems should be hard to simulate is indeed true 
for classical systems, where the frustrated or glassy models 
are the hard ones to describe. For quantum systems, there is 
evidence that the situation should be more complex ETI . 

In this work, we explore a class of models where the in- 
tuition of frustration-free models being easy to solve holds 
true. Building upon work in Ref. |22| and in Ref. |23l . for 
a natural class of two-local Hamiltonians acting on spin-1/2 
particles (simply "spins" henceforth), we show that ground 
states can be reduced to a completely characterized and low- 
dimensional subspace, and then re-constructed by identifying 



the ground state-space of each interaction of the Hamiltonian 
term-by-term. Specifically, the ground space is the image of 
a symmetric subspace under an explicitly constructible, and 
efficiently contractible, tensor network. It follows that the 
ground states satisfy an area law, and hence contain little en- 
tanglement in the above sense. This generalizes recent results 
regarding the existence of states which have little entangle- 
ment in the ground-state manifolds of such Hamiltonians ll24l . 
We discuss how to efficiently simulate the ground state mani- 
fold, and suggest how this could be used to simulate "almost" 
frustration-free quantum lattice models. 



II. PRELIMINARIES 

A. Frustration-free Hamiltonians and area laws 

We consider spin-1/2 Hamiltonians on a lattice. The lattice 
is described by some graph, the vertex set of which we de- 
note by V. Naturally, the Hamiltonian will be local, or more 
specifically include only nearest-neighbor interaction terms. 
We represent the Hamiltonian as 

{a,b}CV 

for some terms h a j, acting on pairs of spins in the lattice de- 
scribed by V . By rescaling, we may without loss of generality 
require that the ground state energy of each interaction term 
h a j, is zero. We wish to describe properties of the ground 
state manifold M of such Hamiltonians, given the list of the 
individual two-spin terms h a j, as input. An important class of 
Hamiltonians are those which are frustration-free (or unfrus- 
trated), for which each ground state vector |$) 6 M is also a 
ground state of the individual coupling terms: i.e. for which 

K, b |$) = (2) 

holds for all h a j, and all |$) € M. The actual ground state p 
is the maximally mixed state over M, and so is mixed unless 
the ground state manifold M is non-degenerate. 
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Our main results pertain to frustration-free spin Hamiltoni- 
ans as above, with the further constraint that each term h a .b 
has at least one entangled excited state 



we 



In |Section VI 

show that the ground states p of such Hamiltonians satisfy an 
area law ifTl [41 4T71 : that is, for a contiguous region of spins 
A C V, the entanglement of formation. If the ground state 
is non-degenerate and hence pure, the entanglement of for- 
mation is nothing but the usual entanglement entropy, of the 
ground state satisfies 



E F {p) < C\dA\ 



(3) 



for C > constant, where \dA\ is the "boundary area" (i.e. 
the number of edges in the interaction graph of H, which are 
incident to both A and V \ A): see |Fig. f| The entanglement 
of formation is the largest asymptotically continuous entan- 
glement monotone, so this also implies an entanglement area 
law for e.g. the distillable entanglement. Thus, ground states 
of frustration-free Hamiltonians contain little entanglement. 
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FIG. 1. The entanglement scales at most as the boundary area \ dA\ 
of a distinguished region A in some lattice, here a cubic lattice, with 
tighter bounds for special cases. 

This follows as a consequence of the fact that the ground- 
state manifold is the image of the symmetric subspace on n 
spins (for some n bounded above by the number of spins of the 
Hamiltonian) under an efficiently simulatable, and explicitly 
constructible, tree-tensor network. We show this in |Section V| 
in which we demonstrate how this allows us to efficiently sim- 
ulate the ground space of the spin models we consider. For the 
cases where we have a tree tensor network with a single top 
root, the problem considered here may be viewed as exactly 
the converse problem to the one discussed in Ref. [25 1. 

Our result of an analytical area law complements results on 
area laws in harmonic bosonic systems J7| [TO] [TT], fermionic 
lfT3l4T5l on cubic lattice and general gapped models in one- 
dimensional quantum chains [16]. For a comprehensive re- 
view on area laws — and on implications on the simulatability 
of quantum many -body systems — see Ref. (TJ. 



B. Quantum 2-SAT problem 

The arguments behind our analysis builds upon and extends 
the ideas of Ref. [22 1, which defined the problem of quantum 
satisfiability, and presented Bravyi's algorithm for QUANTUM 
2-SAT. We describe here the connection between this problem 
and frustration in spin Hamiltonians. 



QUANTUM 2-SAT is the quantum analogue of the "clas- 
sical" 2-SAT problem on boolean formulae. The latter 
asks when there exists an assignment of boolean variables 
xi,...,x n which simultaneously "satisfy" a collection of 
constraints on pairs of those variables. This problem is ef- 
ficiently solvable 11261 : in contrast, the similar problem 3- 
SAT (in which constraints apply to triples of variables) is NP- 
complete J27). In QUANTUM 2-SAT, individual clauses on 
boolean variables are replaced by projectors p a j, with 



Pl,b 



Va.b 



(4) 



on pairs of spins: an instance of QUANTUM 2-SAT is satisfi- 
able if there is a vector which is a zero eigenvector of each 
projector simultaneously. 

For an instance of QUANTUM 2-SAT, determining whether 
there exist such simultaneous zero eigenvectors is equivalent 
to determining whether the Hamiltonian obtained by summing 
the projectors is unfrustrated. Conversely, the problem of de- 
termining when a 2-local spin Hamiltonian H is frustration- 
free may be reduced to QUANTUM 2-SAT, by rescaling the 
terms of the Hamiltonian H so that each term h a ,b has a min- 
imal eigenvalue of 0, and replacing each rescaled term h a b 
with the projector p a (, onto img(/i a .&)- By construction, such 
a substituation does not affect the ground space of the terms. 
Thus, solving QUANTUM 2-SAT is equivalent to determining 
whether a 2-local spin Hamiltonian is frustration-free. 

Recently, random instances of QUANTUM k-SAT with rank- 
1 projectors [28| have been studied for k ^ 2, delineating 
the "boundary" of frustration in fc-local spin Hamiltonians in 
terms of the density of interactions l29]432l . We instead ex- 
tend the findings of Ref. 1221 for k = 2, remarking on implica- 
tions for simulating the ground space manifold in frustration- 
free Hamiltonians. In |Section ITT] we review Bravyi's algo- 
rithm for QUANTUM 2-SAT, in order to demonstrate impor- 
tant features of the reductions involved when they are applied 
to unfrustrated Hamiltonians satisfying natural constraints. 



III. REDUCTION TOOLS FOR FRUSTRATION-FREE 
HAMILTONIANS 

Bravyi's algorithm for QUANTUM 2-SAT [22] efficiently 
demonstrates the satisfiability of an instance of QUANTUM 2- 
SAT by a sequence of reductions of Hamiltonians, yielding a 
homogeneous instance (in which all projectors have rank 1), 
and then verifying the satisfiability of these instances. We 
may similarly use Bravyi's algorithm to detect frustration in 
2-local spin Hamiltonians, and consider the features of these 
Hamiltonian reductions when applied to particular classes of 
frustration-free Hamiltonians. 

Throughout the following, we admit representations of 
Hamiltonians H' with non-zero single spin terms h a , 



H' = J2 <b + ]T h a 

{a,b}CV a£V 



(5) 



and again describe H' as unfrustrated if there exists a joint 
ground state with eigenvalue zero of all terms (including the 
single-spin terms h a ). 



A. Reductions by isometries 

Condensing the analysis of Ref. |22|, we consider a re- 
duction for 2-local Hamiltonians H to Hamiltonians H' on 
fewer spins, provided that H contains only positive semidef- 
inite terms which have non-trivial kernels. Throughout, we 
denote C 2 by V.2- 



1. Two-spin isometric contractions 



Consider a Hamiltonian term h u v of rank 2 or 3. If H 



of di- 



is frustration-free, h u v fixes a subspace of H 
mension at most 2, over which the reduced state p„ „ of a 
state vector |$) £ kei-(H) must be a mixture. We describe 
this reduced state by an encoding of one spin into two. Let 
{IV'o) j l^i) ; 1^2) j IV'a)} be an orthonormal basis for H2 ® %2 
such that 



ker(ft u „) C span {|^ ) , 1^)} . 



(6) 



Define an isometry U u -.u,v '■ 



■U. 



${u,v} 



such that 



E 

cG{0,l} 



ce{o,i} 



(7) 



This is an isometric reduction, similar to those in a tree tensor 
network or aMERA ansatz [33 J . By construction, the support 
of the reduced state p UjV lies in img(f7 u:mI ). We may then 
define a Hamiltonian 



H — U^. UV HU U:U _- 



(8) 



on the subsystem V' = V \ {v}, where the spin v is essen- 
tially deleted; any state vector |$) £ ker(H) then has the form 



\$) = U u:uv \&), \&)£ker{H') 
We may express H ' as a sum of terms 



(9) 



h a ,b ~ U u , u v h at t)U u . MtV , h a 



UU v h' a Uu:u,v. (10) 



(In the case that h is of rank 3, these will include a non- 
zero single spin operator h' u v which acts on u alone.) If H 
contains non-zero terms h a>u and h a . v , we obtain two terms 
^a.u — ^u-.u.v^o-^Uu-.u^v £Uk1 h av — U^.^ yha^Uu-.u^ in the 
Hamiltonian H, which both act on the spins u and a. We sum 
these to obtain a combined term 



h' 



(ID 



in the reduced Hamiltonian, which may be of higher rank than 
either h' a u or h' a v . (We similarly accumulate any single-spin 
contributions h' u and h' v which may arise from single-spin 
terms h u and h v in H.) Fig. [2] illustrates the effect of mul- 
tiple reductions on the interaction graph of the Hamiltonian. 

For example: Consider an anti-ferromagnetic four-spin 
Hamiltonian, with interactions 



2 x u x ' 2 y V 



(12) 





















X 







FIG. 2. Illustration of transformations of the interaction graph of a 
Hamiltonian H by two-spin isometries. Darker, thicker edges repre- 
sent interaction terms h a ,b of rank 2 or 3, which may be eliminated 
by contraction (corresponding to a reduction H h-> U^. ab HU a:a b). 
In the case of an interaction term with rank(/i r ,u) = 3 as illustrated 
above, a single spin operator (represented above by a loop) is pro- 
duced on the contracted vertex. 



on a four-spin cycle with interacting pairs (1, 2), (2, 3), (3, 4), 
and (4, 1). Here we denote 



1 

1 

-i 

1 



|1)(0| + |0)(1|, (13a) 
: »|1){Q|-»|0)(1|. (13b) 



Rescaling the interactions to have ground energy zero (and 
taking these for the hj^ instead) gives us 



h jjh = |0,0)(0,0U + |l,l)(l,ll, fe 



(14) 



The kernel of this operator is clearly spanned by |0, 1) and 
1 1 , 0). We may consider the effect of contracting spins 3 and 
4 into a renormalized spin, using the isometry 

i? 3:3 ,4 = ( |0) 3 ®|l>4)<0| 3 + ( |1) 8 ® |0)4) <l| a • (15) 

By construction, we have i?g. 3 4/13,4^3:3.4 = 0; and as {1, 2} 
is disjoint from {3,4}, we have h! x 2 = R 3 . 3 4/11,2-^3:3,4 = 
hi 2- We compute the renormalized terms h' 2 3 and h' 4 x as 

h'2,3 = 4:3,4 ( 10, 0)<0,0| 2)3 ® 1 4 + |1, 1)(1, 1| 2)3 ® U)R 3 : 3 A 

= |0,0)(0,0| 23 + |1,1)(1,1| 2! 3 , (16a) 

K,l = 4:3,4 ( J 3 ® |0,0)(0, 0| 4)1 + Is ® |1, 1)(1, 1 ;., ) #3:3,4 

= |1,0><1,0| 3 ):L + [0, 1><0, 1| 3 x ; (16b) 
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up to rescaling, the resulting renormalized Hamiltonian is then 

H> = 4:3,4^3:3,4 ~ 

- W+^W]), (17) 



with a ferromagnetic coupling between the site 1 and the 
renormalized site 3. Arbitrary rank-2 or rank-3 interactions 
may be contracted similarly. 



2. Single-spin deletions 

For a 2-local Hamiltonian H containing non-zero single- 
spin operators h v (e.g., such as may arise from the preceding 
reduction), any state vector |$) G ker(iJ) must be factoriz- 
able into a single-spin pure state vector £ ker(/i„) acting 
on v, and some state of the remaining spins. If the operator h v 
has full rank, it follows = 0, so that H has trivial kernel. 
Otherwise, \tp) may be taken to be a unit vector spanning the 
kernel of h v , and we may form a Hamiltonian 



H 1 



l)H(\^) v ®l) 



(18) 



on the subsystem V = V \ {v}, consisting of a sum of terms 
h' a = h a acting on individual spins o/d, and terms 



Va,b=({ll>\ v ®l)ha,b(\ll>)v®l), 



(19) 



acting on pairs of spins {a, b} C V. In the latter case, if 
v e {a, b}, then h' b will be an operator acting on a single 
spin; otherwise, we have h' ab = h at j,. Again, in the case 
of single-spin terms h' a b acting on a spin a, if there was a 
term h a present in H, we may accumulate these into a term 

K = K + h' ab . 

We may again describe H' using an isometric reduction in 
this case: if we define P v = \ip v )<E)ly, then P v is an isometry 
whose image contains any state vector |$) 6 ker(-ff). We may 
then rewrite Eq. (18) as 



H 1 = PlHP v . 



(20) 



expressing it in a form more explicitly similar to Eq. (8) 



3. Remarks on these reductions 

The reductions above correspond to the reductions pre- 
sented in Ref. (22) for instances of QUANTUM 2-SAT which 
contain two-spin operators of rank greater than 1. This allows 
us to reduce QUANTUM 2-SAT to the special case of "homo- 
geneous" instances (in which all terms have rank 1). 

The key feature of both reductions above is that the kernels 
of H and H' are related by isometries, and so have the same 
dimension. If the Hamiltonian H' has any terms of full rank 
(acting on either one or two spins), it follows that the Hamil- 
tonian H' has trivial kernel; then the same holds for H as 



well. If we do not encounter any full-rank terms, each reduc- 
tion produces a Hamiltonian acting on one fewer spins, even- 
tually yielding a "homogeneous" Hamiltonian (extending the 
terminology of Ref. l22l to Hamiltonians in general, including 
those with single-spin terms of rank 1). 

The choice of the reduction at each stage does not matter, in 
the following sense. As long as we have a Hamiltonian con- 
taining two-spin terms of rank at least 2, and which does not 
contain full-rank terms, we may extend any sequence of re- 
ductions to one which terminates with a Hamiltonian H which 
is either homogeneous or contains a full-rank term. In the 
latter case, the original Hamiltonian has trivial kernel, and is 
therefore frustrated; otherwise, we obtain a "homogeneous" 
Hamiltonian whose kernel may be mapped to that of the orig- 
inal Hamiltonain H by a sequence of known isometries. If 
we can solve the homogeneous case, we may then choose the 
reductions according to whichever criteria are convenient. 

Note that this reduction process, from an input Hamiltonian 
H to a homogeneous Hamiltonian, amounts to a tree tensor 
network of isometries (albeit applied to a vector subspace): 
from a temporal top layer defined by a Hamiltonian containing 
only terms of rank 1, one constructs the ground space of the 
full Hamiltonian H by sequential applications of isometries 
with a simple topology. We develop this observation further, 
and remark on implications for simulating the ground space of 
H, 



in 



Section V 



Note also that for a one-dimensional quan- 
tum chain and a sequential contraction, this construction gives 
rise to a sequential preparation of a quantum state and hence 
to a matrix-product state of small bond dimension. 



B. The homogeneous case 

Given a homogeneous Hamiltonian H' (containing only 
terms of rank 1) acting on some system V', consider a col- 



lection of vectors |/3 a ,t>) £ H 



«>{a,b} 



such that 



h a ,b — \0a,b)(0a,b\ 



(21) 



We interpret each two-spin Hamiltonian term h' a b as a con- 
straint on the corresponding two-spin marginal state p ab of 
a state |$) € ker(if), and attempt to obtain additional con- 
straints on pairs of spins a, b € V by combinations of the 
constraints which are already known. Ref. 11221 shows that if 
|$) lies in the kernel of ((3 a ,b\ an d (Pb c\ acting on the corre- 
sponding spins, it also lies in the kernel of the functional 



(22) 



acting on the spins a and c, where l^ - ) oc |0) |1)— 11) |0) is the 
two-spin antisymmetric state vector. We call such a constraint 
(fi'a c \ an "induced" constraint, and use the term induction of 
constraints to refer t o the oper ation on (f3 a b \ and (j3 b c \ which 



Eq. (22) up to a scalar factor. 



gives rise to J in 

For each induced constraint ((3' uv I on a pair of spins u, v, 
we may add a term 



1/3' )(0' I 



(23) 
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to the Hamiltonian H' , obtaining a Hamiltonian H which (by 
construction) has the same kernel as H' . If H' already con- 
tains a term h u v which is not colinear to the induced term 
h UtV , these may be accumulated into a term h u>v whose rank 
is at least 2, and one may apply a two-spin contraction as 
described in |Section 111 A| Otherwise, we may induce fur- 
ther constraints from the terms of H, until we obtain a com- 
plete homogeneous Hamiltonian H c : a Hamiltonian in which 
the two-spin constraints (fi u ,v\ are closed under constraint- 
induction. 

By inducing constraints on pairs of spins, possibly perform- 
ing two-spin contractions as in Section HI A| when we obtain 



terms of rank 2 or more, we may efficiently obtain a com- 
plete homogeneous Hamiltonian H c from a frustration-free, 
2-local Hamiltonian H. Furthermore, Ref. Il22l shows that a 
complete homogeneous Hamiltonian H c acting on at least one 
spin (and which lacks single-spin operators l34l ) has a ground 
space which contains product states. Thus, for homogeneous 
Hamiltonians H, we may either efficiently determine that it is 
frustrated, or efficiently obtain a Hamiltonian which is closed 
under constraint-induction. In the latter case, we may con- 
struct product states in the kernel of H by selecting states for 
each spin consistent with the two-spin constraints 



IV. UNFRUSTRATED NATURAL HAMILTONIANS 

We now present results concerning the ground state man- 
ifold of a "physical" class of 2-local spin Hamiltonians. We 
will say that a Hamiltonian H is natural if it is 2-local, con- 
tains no isolated subsystems, and each term h a j, (acting on 
H2 <8> H2) has at least one entangled excited state (i.e., there 
exists an entangled state orthogonal to the ground state man- 
ifold of h a j,). Without loss of generality, we may further 
require that the ground energy of each term in H is zero. 
This is a natural assumption that typical physical interactions 
will satisfy: for instance, ferromagnetic or anti-ferromagnetic 
Ising interactions (which have excited eigenstates |0)|1) — 
|1) |0) and |0) |0) + |1) |1) respectively), ferromagnetic or anti- 
ferromagnetic XXX models (which also have those respective 
eigenstates), or indeed any interaction which is inequivalent to 
either diag(0, 0, E a , Ef,) or diag(0, E a , 0, E^) up to rescaling 
and a choice of basis for each spin. 

Using the reductions of |Section III A| we show strong 
bounds on the dimension of the ground space of an unfrus- 
trated natural Hamiltonian on spins. This will allow us in 
|Section V] to describe a scheme for efficiently simulating the 
ground space of frustration-free natural Hamiltonians, and 
in 



Section VI to demonstrate that the ground states of such 



Hamiltonians satisfy an entanglement area law. 



Ground-spaces of unfrustrated, natural, homogeneous 
Hamiltonians 



We now present an extension of the analysis of Ref. 
for homogeneous and complete Hamiltonians H c (acting on 
a set V c of spins), to examine the ground-state manifold of 



H c in the case that H c is also natural. We show, using tech- 
niques similar to those used in Ref. 11301 Section III A], that 
the ground space of such a Hamiltonian is equivalent to the 



symmetric subspace Symrr^'Hf ) C Hf c , up to some ef- 
ficiently constmctible choice of invertible operations on each 
spin. As we may reduce more general Hamiltonians, i.e. hav- 
ing terms of rank 2 or 3 (extending beyond those Hamiltoni- 
ans considered in Ref. 11301 ) to homogeneous natural Hamil- 



tonians via the reductions of Section III A these results yield 
important consequences for natural frustration-free Hamilto- 
nians in general. 

Consider a Hamiltonian H c acting on V c , where H c has 
no single-spin terms. Because the two-spin constraints de- 
scribed by the terms of H c are closed under the induction of 



constraints (as described by Eq. (22) 1, and as there are no iso- 
lated subsystems, it is easy to show that every pair of spins 
is acted on by a non-zero term in H c . For such a Hamilto- 
nian, the excited states \0 a ,b) for the terms /i a (, in H c are en- 
tangled states. We may then construct a family of operators 
{L v } veVc Q GL(2) such that 



oc (* \ uv {L u ®L v ) 



(24) 



for each pair of spins u, v, where is again the two-spin 
antisymmetric state vector. For instance, one may fix L a = 1 
for an arbitrarily chosen spin a G V c , and determine linear 



operators L v satisfying Eq. (24) for each v E V c and opera- 
tor ((3 a ,v\- Any such choice of operators {L v } v£V satisfies 
Eq. (24)| for all u, v, which follows from the closure of the 
constraints (/3 U>V \ under induction: 

(Pu,v\ « (<A,«|®<A,,,,|)(1<8>|¥-)<8>1) 

((*-| <8>(#-|)(£ u <g> |*")<E> L v ) 
(*-| (L U ®L„). (25) 



oc 
oc 



We define scalars A Uj „ 7^ such that (j3 uv \ 
^u,v{^~\ (L u ® L v j for each u,v € V c , and let T 
1 : we then have 



(®„ 6 v. L ") 



({P U , V \®1)T 



where we define by 



Tjr 



l: 



(26a) 



(26b) 



As the operators T^j have full rank, the operator ((j3 UiV \ 



1)T then has the same kernel as ( , 3>~ 



, which is the +1- 



eigenspace of the SWAP operator acting on u and v. It follows 
that the kernel of the Hamiltonian 

T^H C T = Tt(|/3 u ,„)(/3 u ,„|(g)l)r 

u,vGV c 



is the symmetric subspace Symm('Hf 1 
the result of J30] Section III A]. 



: ) ; this corresponds to 
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We remark on some important properties of Symm('Hf mc ), 
where n c = \V C \. This subspace is spanned by uniform super- 
positions \Wk) of the standard basis states having Hamming 
weight ^ k ^ n c , 



\W k ) cx J2 



(28) 



x<£{0,l} r 



thus dim Symm("H2 



n c + l. This subspace may also be 



spanned by product state vectors \ao 



K) 8 "' for 



any set of n c +l pairwise independent state vectors \af) 6 H.2- entangled excited states. As h' a 



for orthonormal two-spin state vectors l^o) , l^i) whose span 
contains ker(/i u „); we require that \ipi) be entangled, which 
ensures that U u:UtV is not a product operator. Consider the 
terms h' a b — U\. u v h a ^U u:UtV which occur in the Hamilto- 

nian H' = Ul. u v HU u:UtV . For any two-spin operator h v ^ a 
acting on v and some other spin a, the fact that h v<a has an en- 
tangled excited state implies in particular that it is not a prod- 
uct operator. Thus, h' v a is a product operator only if it has full 
rank. If H is frustration-free, h' v a cannot have full rank; then 
h' v a is not a product operator, and in particular it will have 

h a b when a, b ^ {u, v}, 



Thus, any natural Hamiltonian H c which is also complete and 
homogenous has a ground space of dimension n c + 1, and can 
be spanned by a family of classically efficiently simulatable 



state vectors 



(L v \a>j)), 



(29) 



v£V c 



for some choice of pair-wise independent single-spin state 
vectors |ao) , ■ ■ ■ , low) we use this fact in Section V 



Note that if even this efficient method should be too compu- 
tationally costly for very large systems, one can also Monte- 
Carlo sample from the ground state manifold in this way. 



B. Preservation of natural Hamiltonians under reductions 



A key feature of natural Hamiltonians (defined on page 5 1 



is that the class of frustration-free natural Hamiltonians on 
spins is preserved by the two-spin contractions described in 
|Eq. (8)| This implies that the reductions of |Section III A| map 
the ground-state manifold of an unfrustrated natural Hamilto- 
nian H provided as input to that of a complete, homogeneous, 
natural Hamiltonian; we may then apply the results of the pre- 
ceding section to describe the ground-state manifold of H . 



Consider an isometry U u 



n 



«>{«} 



h: 



de- 



rived from a two-spin Hamiltonian term h u , v as described in 
Section III A 1 We may show that for any term h a u in H, the 
corresponding term h' a u — £7j. u v h atU U u:U!V in the reduced 
Hamiltonian H' = U^. u V HU U:U ^ V has an entangled excited 
state if the same holds for h a ^ u . We require the following two 



lemmas, whose proofs we defer to Appendix A 



Lemma 1 (Product states). For two-spin state vectors \ip) and 
\4>), we have {{tp\ ® t) (l (g) \(j>)) = only if both \ip) and \<j>) 
are product states. 

Lemma 2 (Product operators). Let U : % 2 — > H2 <& be 
an isometry which is not a product operator. Let r] ^ be an 
operator on two spins, and 7/ = (W ®l2)(~i-2®v)(U <8>J-2)- 
// 77' is not of full rank, then rj is a product operator if and 
only if r] is a product operator. 



it follows that H' is a natural Hamiltonian; and as H' has a 
kernel of the same dimension as H, it is frustration-free as 
well. 

We may strengthen this result, to show that if H is natu- 
ral and frustration-free, and also contains no two-spin terms 
of rank 1, then the same is true of H' = U^. u V HU U:U:V as 
well. For any two-spin term h v a acting on v with rank at least 
2, consider states \(po) , \tpx) € img(/i„ ia ) such that \ipo) — 
\a) is a product state and \ipi) is entangled; Any subspace 
of %2 <8> H2 °f dimension at least 2, such as img(/i u „), con- 
tains a product state vector \>po); the existence of \<pi) is guar- 
anteed by the definition of a natural Hamiltonian (compare 
also Ref. l37l ). and choose real parameters Ao, Ai > such 
that 



h v , a ~ A l^oK^ol - Ai \ipi)(ipi\ ^ 0. 



(31) 



Let rjk = \fk)(<Pk\ for k € {0, 1}, and consider the images i]' k 
under contraction by U u - UtV : 



Vk 



= EdiXV'ji ® 1) (1 ® {m(£\ 8 1) 



J2\j)(t\®M 3 , k Ml 



(32) 



where we define Mj^ = ( (ipj \ ® l) (l (8> \ifik) ) • By Lemma 1 



we have Mj.k = only if both \ipj) and \ipk) are product 
operators; this implies that the operators M± k m particular 
are non-zero, so that rj' k ^ for any k. Note that 



h'v, a ~ ^u-.u.v^v^Uu-.u^ 
:u,v 

= A ?7 + Air/i ; 



AxC/t 



(33) 



because H' has a non-trivial kernel, h' v a has rank at most 3, 
in which case neither operator rj' k has full rank. BylLemma 2 



We show that frustration-free natural Hamiltonians are pre- rank at least 2; by Eq. (33) 



7/4 is not a product operator; as 77 is a product operator, these 
operators are linearly independent. Then, XqT]' + Ai?^ has 

the same is true of h 



served by the reductions of Section III A 1 as follows. Let H 



be a natural 2-local Hamiltonian, and h u , v be a two-spin term 
in H. Define 



If all of 

the terms in H have rank 2 or higher, the same then holds for 
H' as well. 



U u:u<v = |7/> )(0| + foh)<l| 



(30) 



Thus, if we apply the reductions of Section III A to an initial 
Hamiltonian which is both natural and frustration-free, the re- 
sulting Hamiltonians will also be natural and frustration free. 
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Furthermore, if H contains no terms of rank 1, then neither 
will the reduced Hamiltonians. Because the process of induc- 



ing constraints described in Section MB also preserves the 
property of each term h a j, having entangled excited states, 
these invariants ensure that initial Hamiltonians with these 
properties (natural and frustation-free, and possibly contain- 
ing no terms of rank 1) may be reduced to homogeneous and 
complete Hamiltonians which have these same properties. We 



may then apply the results of Section IV A to these reduced 
Hamiltonians. 



P 



t 



r t 


T T 






1 


1 









FIG. 3. A simple tree tensor network. 



V. SIMULATING GROUND SPACES OF 
FRUSTRATION-FREE NATURAL HAMILTONIANS 



Building on the results of Section IV we now show how 



the reductions of Section III A may be used to obtain a proce 



dure for simulating states from the ground-state manifold of 
frustration-free natural Hamiltonians H. 



(neither directed nor undirected): the output indices of each 
tensor represent spins whose state depends on only a single 
spin at the input. Put another way: any spin v which is intro- 
duced by an isometry U u - U . v may be considered a "daughter 
spin" of a unique parent U, which imposes a tree-like hierar- 
chy on the tensor network T, as illustrated in |Fig. 3| Strictly 
speaking, the quantum circuit or tree tensor network will have 
the structure of a forest graph, which is a graph which may 
have more than one connected component, each of which are 
trees. 

The roots of each tree are spins u which are either prepared 
by an isometry P u derived from the removal of single-spin 
terms, or which correspond to free indices at the input of the 
tensor network T. 

In the case that H is non-degenerate, the resulting tensor 
network will (by that fact) simply be a tree-tensor network 
with no free input indices. Conversely, if the input Hamilto- 
nian H is degenerate, there will necessarily be free input in- 
dices, representing a domain consisting of a state space of di- 
mension at least 2. In the latter case, the tensor network T will 
yield ground states of the original unfrustrated Hamiltonian H 
if and only if it operates on a state \cp) E kei(H c ) at the input, 
where H c is the complete homogeneous instance obtained by 
the Hamiltonian reductions. Thus, if one may efficiently sim- 
ulate states from the ground space of such a Hamiltonian, we 
may apply the network T to simulate the ground space of the 
original Hamiltonian H. 



Efficiently simulating ground spaces of unfrustrated 
natural Hamiltonians 



A. Tree tensor networks and matrix-product states 



As we noted on page 4 implicit in the reductions of Sec- 



|tion III A| is that the isometric reduction from general Hamilto- 
nians to homogeneous instances has the form of a tree-tensor 
network. Thus, simulating the ground state manifold of any 
unfrustrated 2-local Hamiltonian on spins may be reduced 
to that of a complete homogeneous Hamiltonian acting on a 
smaller system. In this section, we sketch this reduction. 
For any unfrustrated 2-local Hamiltonian H on n spins, we 



may apply two-spin reductions as described in Section III A 



until we obtain a homogeneous instance without single-spin 
terms. We then attempt to induce additional constraints via 
|Eq. (22)| and apply further two-spin reductions if we obtain 
terms of rank 2 or 3. If H is frustration-free, this process will 
ultimately terminate in a complete homogeneous Hamiltonian 
H c on a subset V c C V. 

Consider the tensor network T which performs the com- 
plete reduction as above; we describe T in reverse order, as 
introducing new spins to represent a unitary embedding of 
ker(i? c ) into kcr(H). The various spin contraction isometries 
U u -.u,v as in Eq. (7)| each have a single input-index and two 
output indices; the spin deletion isometries P v have no input 
indices at all. These are applied sequentially, giving rise to an 
acyclic directed network. As the in-degree of each tensor is 
at most 1, it follows that the network contains no cycles at all 



Tensor networks T with free input indices, and with a tree- 
like structure such as described above, can be efficiently sim- 
ulated over inputs with low Schmidt measure [38], as follows. 

For any observable £1 acting on m spins, one may eval- 
uate (f2)„ for the maximally mixed state over the ground- 
state manifold of H by computing the expectation (T^DT') H 
over the ground-state manifold of the complete homogeneous 
Hamiltonian H c obtained as described in |Section III| As the 
tensor T has tree-stucture, the observable also acts on 

at most m spins. If we can obtain an orthonormal basis for 
ker(iJ c ) which may be succinctly described in terms of prod- 
uct states, we may evaluate expectation values of T^ilT with 
respect to m-fold products of single spin states. 

As we noted in Section IV A| ker(_ff c ) can be spanned by 
a collection of n c + 1 product vectors (where n c is the num- 
ber of spins on which H c acts). Let \$' ) , , . . . , \& nc ) € 
ker(_ff c ) be a collection of independent product vectors, 



(34) 



We may efficiently compute a projection of il c onto ker(_ff c ) 
by performing a suitable transformation of the matrix 



W(Q C 



E 

j,k=0 



\j)(&j\il D \tf k ){k\ 



(35) 
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as follows. The operator in particular is positive defi- 

nite; we thus have W^(ll) = U AW for some unitary U unitary 
and positive diagonal matrix A. It is not difficult to show that 



j=0 



3=0 



(36) 



for some orfhonormal basis |$q) , ■ ■ ■ > \ ®n c ) of ker(H c ), by 
taking the product of the above operator with its adjoint. Thus, 
the restriction of fl c to ker(H c ) with respect to the basis of 
states may be computed as 



n = a- x / 2 Ww{q. c )ua- x 



/2 



(37) 



By considering operators fl c = T^QT, this allows us to com- 
pute the restriction of operators to the ground-space of Hjj: 
seef 



Fig- 4 



A. Area law for frustration-free natural Hamiltonians 

We consider first the case where A is a contiguous subsys- 
tem (i.e. for which there is a path between any pair of spins 
in A, following edges in the interaction graph of the Hamilto- 
nian H), and subsequently generalize the observation in this 
case to arbitrary subsystems A. We first decompose 



H = H/ 



H B + H_ 



(38) 



for B = V \ A, and where Ha and H B contain all terms 
internal to A and B, respectively. We then apply the reduc- 
tions of Section III A to the subsystem A. That is, we per- 



form two-spin contractions as described in Section III A 1 for 
any two-spin terms in Ha of rank 2 or 3, and perform spin- 



deletions as described in Section III A 2 for any single-spin 



terms in Ha- Performing the constraint-induction process of 
Eq. (22) — again on the terms acting on A alone — and then 









a Hamiltonian H of the form 








H = H A + H B + H Ab ' 


Tree tensor network 



JUL 





Tree tensor network 




1 1 






Symm(7C=) 





FIG. 4. Schematic diagram showing the isometric decomposition for 
efficient simulating the ground state manifold. A local Hamiltonian 
term supported on two sites is marked yellow. 

We may thus efficiently estimate such observables with re- 
spect to ground-states of Hjj\ for constant m, the required 
inner products may be calculated as (sums of) scalar products 
of at most n c inner products over vector spaces of bounded 
dimension. To evaltuate the value of ft with respect to ground 
states of the input Hamiltonian H, it suffices to analyse the 
polynomially-sized operator Ct representing the action of ft 
on the ground-state manifold. 



VI. ENTANGLEMENT BOUNDS FOR GROUND STATES 
OF FRUSTRATION-FREE NATURAL HAMILTONIANS 



The fact that the reductions of |Section III A preserve the 
class of natural Hamiltonians (as defined on page 5 i allow us 
also to make more global statements about ground states for 
frustration-free Hamiltonians, again by reduction to the com- 
plete homogeneous case described in |Section IV A| In this 
section, considering frustration-free natural Hamiltonians H, 
we demonstrate an area law for the entanglement possible in a 
ground state of H between any subsystem A C V and its en- 
vironment V \ A. We also consider some very general cases 
in which still stronger upper bounds on the entanglement may 
be obtained. 



(39) 



where A C A is the set of spins remaining after the reduc- 
tion process, and where H A is a homogeneous and complete 
Hamiltonian. The Hamiltonians Hs and Hz B contain the 
terms derived from Ha and Ha.b respectively by the reduc- 
tion process. In other words, we perform a complete tree ten- 
sor reduction on the subsystem A, until we obtain a Hamil- 
tonian whose restriction to A is homogeneous and complete. 
As dimker(if) = dimker(iJ) > 0, we have ker(ii^) > 
and ker(H B ) > as well; in particular, ker(H) C ker(iTj), 
so that Hx is frustation-free. Because ker(if) C ker(H B ) as 
well, we then have 



ker(H) C ker(H A ) <g> ker(H B ) 



(40) 



taking the restrictions of Hz and H B to their respective sub 
systems A and B. Let h = \A\: as H% is also homoge 
neous and complete, it has kernel of dimension h + 1 by Sec 
Ition IV Al 

In the case where A contains multiple components 
A\, A-2, . . . , Ak with respect to the interaction graph of H 
(where each Aj is disconnected from the others but connected 
internally), we may perform the Hamiltonian reductions of 



Section III A to each component independently. We may f ur- 
ther decompose the Hamiltonian Hi obtained in ] 



Eq. (39) 



H A = 



Hi 



where Aj = A n Aj. (41) 



As H 7 is unfrustrated, each of the sub-Hamiltonians H i is 
unfrustrated as well, in which case we may write 



ker(i^) C ker(H Ai ] 



>ter(H Ak ), (42) 



similarly to Eq. (40) 



Then, the dimension of ker(H A ) is 
bounded by the product of dim ker(H A ) = hj + 1 for each 
subsystem, where hj — \Aj\. Let aj be the number of spins 
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in Aj which are adjacent in the spin lattice to spins in B, and B. Stronger entanglement bounds for contiguous subsystems 

let o.j ^ ay be the number of such spins in Aj. For any H 
with nearest neighbor interactions on a lattice in finitely many 
dimensions (in the graph theoretical sense), there exist scalars 
c, K > such that 



\j > Kn C j 



(43) 



for each subsystem Aj. We may then bound on the dimension 
of ker(H^) in terms of these "boundary spins" as 

k 

log(dimker(ff y |)) < ^ log(dim ker(jfj . )) 



3 = 1 



3 = 1 



(44) 



- ak is the number of spins in A adjacent 



where a = a% + ■ 
to elements of B. 

From this bound on dimker(i?^), it follows that any pure 
state |$) in the ground space of H has Schmidt measure l38l 
at most a/K across the partition A + B, and so can support 
at most this many e-bits of entanglement between A and B. 
Because any state vector |\I>) £ ker(H) can be obtained from 
some vector |$) £ ker(H) by a network of isometries act- 
ing only on spins in A, it follows that any state vector |^) 
in the ground-state manifold of H contains at most a/K e- 
bits of entanglement between A and B. Note the similarity to 
ground states of low Schmidt rank close to factorizing ground 
states in Heisenberg models [36]. In case the ground state p is 
degenerate, each pure state in the spectral decomposition of p 
will have that property. As the entanglement of formation is 
convex (this usually being taken as a necessary property of an 
entanglement monotone), one obtains the bound 



E F (p) s? 



K 



(45) 



for the entanglement of formation between A and B. Finally, 
let dA be the set of edges between A and B. By definition, 
for each edge in dA, there is a spin in A which is adjacent to 
some spin in B; then we have a ^ \dA\, so that 



E F {p) < 



K ' 



(46) 



Thus, the amount of entanglement which can be supported by 
a ground state of H between A and B is governed by an area 
law. We summarize: 

Proposition 1 (Area law). Let H be an unfrustrated natural 
Hamiltonian on a lattice V, and denote with p its (possibly 
degenerate) ground state. Then for any subsystem A C V, the 
entanglement of formation of p with respect to A and V x A 
satisfies an area law, i.e. there exists a constant C > of the 
lattice model such that 



E F (p)^C\dA\. 



The above analysis imposes no additional constraints, be- 
yond the requirement that H be natural and frustration-free. 
We may obtain still stronger bounds — by the logarithm of 
the system size, or even by a constant — on the entanglement 
between A and its environment B — V \ A, under fairly gen- 
eral conditions on the subsystem A when it is a contiguous 
subsystem. 



1. Contiguous subsystems in general 

Implicit in the analysis of the previous subsection is a 
stronger entanglement bound for contiguous subsystems in 
general: we observe that if A consists of a single component, 
we have 



dimker(£f|) = n + 1 



(48) 



for n — | A | by the analysis of Section IV A (where H% is the 
reduced Hamiltonian acting on the subsystem A described in 



Eq. (39) i. By a similar analysis, if a is the number of spins in 
A adjacent to at least one spin in B, we may use Eq. (43) to 
obtain 



log(dimker(_ff^)) = log(n + 1) ^ 
Asa ^ \9A\, we may then obtain: 



log(a/K + 1) 



(49) 



Proposition 2 (Logarithm law for contiguous systems). Let 
H be an unfrustrated natural Hamiltonian on a lattice V, and 
denote with p its (possibly degenerate) ground state. There 
then exists a constant C > of the lattice model such that, 
for any contiguous subsystem A C V, the entanglement of 
formation of p with respect to A and V \ A satisfies 



E F (p)^C\og\dA\. 



(50) 



2. Subsystems acted on by many high-rank Hamiltonian terms 

In the above result, we have neglected the difference in the 
sizes of the subsystem A, and the reduced subsystem A C A. 
The difference in their sizes will be precisely the number of 
isometric reductions performed to obtain A from A. Each iso- 
metric reduction corresponds to either an edge contraction in 
the interaction graph G of the Hamiltonian H, or a vertex dele- 
tion in G, yielding an interaction graph G' for the Hamiltonian 
H' . Two-spin isometries U u:uv represent the reduced state- 
space of the two-spin subsystem {u, v} as the image of a sin- 
gle spin under an isometry: the corresponding reduction may 
thus be represented as a "contraction" of two spins into one, 
as illustrated in Fig. 2 Single-spin terms h u may be repre- 
sented by loops on vertices: isometric reductions arising from 
terms h u>v of rank 3 also yield a loop on the contracted ver- 
tex. Spin-removal reductions P u correspond to the deletion of 
a vertex u with a loop, which removes all edges au incident 



(47) 

to u, possibly replacing them by loops on the neighbors a. 
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This representation of Hamiltonian reductions in terms of 
graphs is underdetermined, in that it is not always possible 
to determine the ranks of the reduced Hamiltonian H' from 
those of the Hamiltonian H prior to contraction. However, 
the correspondence to graph reductions motivates a simple ob- 
servation. Consider a subsystem A, and consider the Hamil- 
tonian Ha together with its interaction graph G a- We may 
"colour" or "rank" the edges of Ga according to whether the 
term corresponding to each edge is rank-1 (which we call 
"light" edges) or has rank 2 or 3 (which we call "heavy" 
edges). The two-spin isometric reductions of Section III A 1 



required to obtain Hi corresponds to contractions of all heavy 
edges in Ga- As such contractions preserve connectivity, this 
implies that the interaction graph Gj corresponding to Hj 
has as many vertices as there are connected components in 
the "heavy subgraph" of Ga- In particular, if the number of 
"heavy" connected components (components connected only 
by heavy edges) is bounded above by some parameter (3, we 
then obtain h = \A\ ^ j3, so that 



log(dimker(iTj)) < log(/3 + 1) 



(51) 



A consequence of this is that if H is a frustration-free nat- 
ural Hamiltonian which contains only terms of rank 2 or 3, 
all edges in Ga will be heavy, so that it consists of a single 
heavy component; we then have ^(dimkei^iJ^)) = 1. In 
this case, there is at most one e-bit of entanglement between 
A and any other, disjoint subsystem in the lattice. 

We may further refine this observation by considering the 
impact of Hamiltonian terms of rank 3. More generally, we 
may consider rank-1 single-spin terms in the reduced Hamil- 
tonians, arising either from rank-1 terms in preceding Hamil- 
tonians, or from performing an isometric reduction on terms 
of rank 3. Such single-spin terms correspond to loops on ver- 
tices in the interaction graphs G' A , of the reduced Hamilto- 
nians H' A i ■ In the case of a frustration-free natural Hamilto- 
nian containing such terms, we may show that the Hamilto- 
nian H is non-degenerate with a ground state consisting es- 
sentially of a product of single-spin states (together with a 
single two-qubit entangled state if the original Hamiltonian 
H contains a term of rank 3). Consider the effect of preferen- 
tially performing single-spin deletions in the process of reduc- 
ing H by isometries: for a natural Hamiltonian, we may eas- 
ily verify that removal of such a vertex u (i.e. performing the 
single-spin removal reduction of |Section III A 2| > will induce 
loops corresponding to single-spin terms on all neighbors of 
u. These spins may then be removed in turn, inducing still 
further loops; by the requirement that the original interaction 
graph G be connected, this ultimately results in the removal 
of every spin on which H acts, each by independent single- 
spin isometries which describe a fixed single-qubit state vector 
\(p u ). As a result, the entire lattice contains no entanglement, 
or at most one e-bit if H contained a single rank-3 term giv- 
ing rise to a "seed" loop; any subsystem A which is not acted 
on by the rank-3 term therefore contains no entanglement, nor 
has any entanglement with its environment. 

The above exhibits the fragility of the condition of 
frustration-freeness: it follows, for instance, that any natu- 
ral Hamiltonian H which contains as many as two terms of 



rank 3 is necessarily frustrated (i.e. does not have a ground 
space characterized by those of its interaction terms). Because 
the same unique ground state must be produced by any re- 
duction, e.g. in which we first perform two-qubit isometries, 
it follows that each two-qubit isometry in such a reduction 
must also map the single-spin states (describing the unique 
ground state of the reduced Hamiltonians) to product states, 
which is of course highly unlikely if instead one considers 
arbitrarily chosen two-qubit isometries and single-product in- 
put states. These observations may be used together with the 
random satisfiability results of Ref. 1 30 1 to suggest that "ex- 
act" frustration-freeness is likely to be rare in physical sys- 
tems; small perturbations are likely to cause frustration. This 
is nothing but a manifestation of a fragility against sponta- 
neous symmetry breaking. However, in Section VIII we sug- 



gest ways in which systems which differ only slightly from 
frustation-free systems may be examined using the techniques 
of lSections"Vl and VII 



VII. DIFFERENT MODELS OF FRUSTATION-FREE 
HAMILTONIANS 

In this section, we consider frustration-free Hamiltonians 
H, but suspend our earlier restriction to natural Hamiltonians 



(as described on page 5 i in order to consider different mod- 



els of Hamiltonians that are of interest. In doing so, we will 
compare the resulting analysis to the case of frustration-free 
natural Hamiltonians in lSection VII 



A. Rank-two terms lacking entangled excited states 

Any Hamiltonian term h a ^m H which has rank 2 and has 
only product states orthogonal to its ground space is of the 
form h a j, — \<p)((fi\ a <8> (or the reverse tensor product), 
where rj is a single-spin operator of full rank and \<p) £ %2- 
This operator has the same kernel as the single-spin operator 
\<p)((fi\ a £?> 1&; therefore, if H is frustration-free, we may per- 
form this substitution without any change to the ground state 
manifold or its properties. As any rank-3 operator has entan- 
gled states orthogonal to its (unique) ground state, we may 
therefore restrict to the case where "non-natural" terms h a ,b 
occuring in H have rank 1, so long as we permit input Hamil- 
tonians with single-spin terms. 



B. Unfrustrated translationally invariant Hamiltonians 

Consider a frustration-free Hamiltonian H in which the 
interaction terms h a j, of each spin a is the same for all of 
its neighbors b. If H is not natural, it follows that h a j, — 
\a)(a\ a £g> |/3)(/3| fc for some states \a) , \(3) € "H 2 ; and by a 
suitable choice of basis on each site, we may without loss of 
generality let \a) = \/3) = |1). 

Consider a ground state vector |$) of the Hamiltonian. For 
each site a, if the state vector of \a) is not given by |1), it 
follows that all of the neighbors of a are in 1 1); and conversely, 
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if all of the neighbors of some site a are in |1), the site a may 
be in an arbitrary single-spin state without contributing to the 
energy of the global state. It follows that the ground-state 
manifold of H consists of all superpositions of product states 
in which all sites are in |1), except for some set of mutually 
non-adjacent sites ACS, whose spins may have arbitrary 
states (including states which are entangled with other sites 
in S). In particular, for bipartite lattices, this includes states 
in which the entire "even" sublattice of sites an even distance 
from the origin are in |1), and the opposite "odd" sublattice 
may have an arbitrary entangled state. 

Thus, if H is isotropic and frustration-free, then without 
loss of generality it is either natural, or contains subspaces in 
which large subsystems of the lattice are essentially uncon- 
strained, and may occupy states with arbitrarily large entan- 
glement content. Consequently, one may expect that unfrus- 
trated translationally-invariant Hamiltonians should have in- 
teraction terms with entangled excited states, i.e. be given by 
natural Hamiltonians. 



C. Unfrustated lattices with randomly located product terms 
and percolation 

Finally, we wish to consider a class of random Hamil- 
tonians which includes non-natural Hamiltonians, and com- 
pare the behaviour of their ground-state manifolds to natural 
Hamiltonians. If one distributes random Hamiltonian terms 
h a _b over nearest-neighbor pairs in an arbitrary lattice, then 
they will almost certainly have an entangled excited state, as 
the highest-energy eigenstate of each term will be a product 
state with probability zero. This remains true even if one con- 
strains each interaction term h a b in the lattice to have ranks 
described by integers r j, € {1,2,3} selected according to 
any distribution, including the case where every term has rank 
1. In order to obtain a random model of non-natural Hamilto- 
nians, we must explicitly designate certain interactions h a b to 
be rank-1 product operators (non-natural terms of higher rank 
being subject to the remarks of |Section VII A| above), and con- 
sider the scaling of the resulting lattice model. 

Consider a d-dimensional rectangular lattice, in which each 
term h a j, has rank-1, and for each term we randomly deter- 
mine whether h a f, is a product term (i.e. satisfies h a & = 



*>« ® for 



some \a) 



£ H2) or an entan- 



gled term (satisfies h a ^ = \^)(^\ for some entangled I7) £ 
H2 (8 The probability that h a j, is entangled is given 

by some fixed p > 0, independently for each edge. Hav- 
ing determined whether h a ,b is entangled or not, we select a 
random rank-1 projector for h a ,b subject to that constraint on 
h a j,. Considering only frustration-free Hamiltonians H con- 
structed under such a model ATI and subsystems A of the lat- 
tice on which H acts, we wish to determine the dimension of 
the ground-state manifold M — ker (Ha)', by a similar anal- 
ysis as in Section VI this will indicate how close the ground- 



states of H come to obeying an entanglement area law. 



As we noted in Section IV A the process of inducing rank-1 
constraints as in Eq. (22) will yield entangled ("natural") con- 
straints from two other entangled constraints. Consider the 



subgraph of the lattice consisting of entangled constraints: it 
follows that any subsystem A C V of the lattice which is con- 
nected only by entangled constraints forms a subsystem for 
which Ha is a natural Hamiltonian, with a kernel of dimen- 
sion at most I A | + 1. Conversely, we may easily show that for 
any product term h a , u = \a)(a\ a ® \P)(P\ U , the constraints 
h a>v induced by h a ,b together with any other constraint h u<v 
will also be a product term (regardless of whether h u <v is a 
product term). Thus, such product terms h a u in the Hamilto- 
nian represent obstacles to the induction of constraints which 
would yield bounds on entanglement: as we noted in Sec 



|tion VII B] above, the prevalence of product terms in a Hamil- 
tonian H allow for the effective decoupling of large subsys- 
tems in the ground-state manifold of H, yielding extremely 
high degeneracy. 

These observations suggest an approach to bounding 
dim(Af) using percolation theory |39| to bound the number 
and size of components connected by entangled edges in a 
large convex subset (for a review on applications of percola- 
tion theory in quantum information, see Ref. {401 ). We may 
consider the worst case scenario in which no additional con- 
straints may be induced between any two subregions A\ , A2 
which are internally connected by entangling terms, but sepa- 
rated by a barrier of product terms which effectively decouple 
the subsystems A\ and A2. If the probability p is above the 
percolation threshhold p c of the lattice, we may apply the fol- 
lowing results: 

Proposition 3 ([39, Theorem 4.2]). Let A be a hypercube 
consisting of n vertices, in a d- dimensional rectangular lat- 
tice with edge-percolation probability p. Then there exists a 
positive real K p eM such that the number of connected com- 
ponents in A grows as n p n, as n — > 00. 

Proposition 4 (|39 Theorem 8.65]). Let C be a finite- size 
connected component containing an arbitrary vertex (e.g. 
the origin) in a d-dimensional rectangular lattice with edge- 
percolation probability p. For p > p c , there exists a positive 
i) p eR such that 



Pr(\C\ = 



exp 



,{d-i)/d 



(52) 



Both k p and exp(— T) p ) in the propositions above are an- 
alytic for p > p c , and thus must converge to in the limit 
p — > 1. We may thus describe an upper bound on the dimen- 
sion of M as follows, for A a large cube containing n ^S> 1 
spins. For p > p c , there is almost surely a unique maximum- 
size component A of the lattice which is connected by entan- 
gled edges: because the percolation probability 6 p is strictly 
positive (by definition) forp > p c , we will have \Aq\ = 6 p n 
on average. Each subsystem Aq, A\, . . . , Ak C A which is 
connected by entangled edges induces a natural Hamiltonian 
Haj which has a kernel of dimension at most \Aj\ + 1: we 
may bound dim(Af) by noting that 



M Cker(i? Ao )(g)ker(iy A J 



>ker(H Ak ), (53) 
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as in Eq. (42) This allows us to obtain the bound 



logdim(M) ^ logdim(ker(if Aj )) 

< E iog(i^i + 1). 

3^0 



(54) 



By Proposition 3 the expected number of components k 



grows like n p n for some n p > as n — > oo; using the prob- 



ability bound on the typical finite component size of Proposi 



|tion 4| as the probability of an indistinguished connected com- 
ponent having size s, we obtain the upper bound 



Exp logdim(M) 



< log(|^ | + l)+^log(|^| + l) 

K p n 

^ \og(e p n + 1) + E 



E 



log(s + 1) 



' exp (?7ps( d - 1 )/ d ) 



log(# p ri + 1) + KpC p n , 



(55) 



where C p is the sum in square brackets (which is small for 
exp(— 7] p ) small). 

As logdim(Af) is also the logarithm of the maximum 
Schmidt rank of any state with respect to the bipartition into 
A and V \ A for the lattice V, the amount of entanglement 
scales with the logarithm of the size of the cube A, with a 
small and tunable linear correction, for p ps 1. In this sense, 
frustration-free Hamiltonians in such a "percolated" product- 
model on rectangular lattices resemble frustration-free natural 
Hamiltonians in the expected case as p — > 1. 



VIII. ALMOST FRUSTRATION-FREE HAMILTONIANS 

The method of efficiently simulating ground state mani- 
folds of frustration-free Hamiltonians can be extended to serve 
as a method to simulate almost-frustration-free Hamiltonians, 
albeit in a non-certified way. Consider a Hamiltonian 

H = H + Ai2i (56) 

for A 6 R playing the role of a small perturbation, where 

(57) 



Hq — E ^ a ' b 

{a,b}QV 



is exactly frustration-free (i.e. in the sense defined in Sec 



tion II A i, and Hi is a small local perturbation. Then, one can 
still efficiently compute 



where M denotes the (in general, degenerate) ground state 
manifold of Ho- Again, we may characterize M as the im- 
age of the low-dimensional subspace Symm(%®" c ) under a 
tree-tensor network, as described in |Section Vj Hi being a 
local Hamiltonian, each term of the infumum above can be 
efficiently computed using a suitable basis of Symm('Hf mc ). 
This is a variational approach that will always provide an up- 
per bound to the true ground state energy. 

In this way, one approximates the ground state manifold 
of an almost frustration-free Hamiltonian with the ground 
state manifold of an exactly frustration-free one. The inter- 
esting aspect here is that one can consider the image of an 
entire large subspace under a tensor network. In practice, one 
would think of a Hamiltonian Hjj near to a realistic one H, 
where one may show that Hjj is frustration-free (which may 
be efficiently verified using the algorithm of Ref. [22 ], as out- 



inf 

|*)eM 



(58) 



lined in Section IIIi, and then approximate the ground state 
of the full Hamiltonian. This approach appears to be partic- 
ularly suitable for slightly frustrated Hamiltonians reminding 
of Shastry-Sutherland type [42] models, with — in a cubic lat- 
tice and a frustration-free Hamiltonian — an additional bond 
along the main diagonal renders the model frustrated. 



IX. SUMMARY 

In this work, we have investigated in great detail a class 
of models whose ground-state manifolds can be completely 
identified: those of physically realistic frustration-free mod- 
els of spin- 1/2 particles on a general lattice. We have seen 
that the entire ground state manifold can be parametrized by 
means of tensor networks applied to symmetric subspaces, by 
essentially undoing a sequence of isometric reductions. We 
also found that any ground state of such a system satisfies an 
area law, and hence contains little entanglement. This is a 
physically meaningful class of physical models — beyond the 
case of free models — for which such an area law behaviour 
can be rigorously proven. It is the hope that the idea of consid- 
ering entire subspaces under tensor networks, and eventually 
looking at the performance when being viewed as a numerical 
method, will give rise to new insights into almost frustration- 
free models. 
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Appendix A: Technical lemmas 

We now supply the proofs of technical lemmata required in 
the preceding sections. 



Proposition (Lemma 1 1. For two-spin state vectors \ifj) and 
we have ( (ip\ <g) 1) (l ® \(j>) ) =0 only if both and \<p) 
are product states. 

Proof. Consider Schmidt decompositions 

|V) = 5> r |e r ) \f r ) ,\<P) = Y1 v ° \Ss) \h.) , (Al) 

where without loss of generality we may require (f r \g r ) ^ 
by an appropriate choice of labels. Then we have 

) 1)(1 ® \4>)) = J2(^s (fr\g s )) \e r )(h s \ , (A2) 



which is only zero if /io^o = Mi^i = 0' which implies fi r = 
and v\_ r = for some value of r £ {0, 1}. ■ 



Proposition ( Lemma 2 1. Let U : H2 — > % ® H2 be an 

isometry which is not a product operator. Let r\ be an 
operator on two spins, and rf = (U^ ® I2XI2 ®tf)(U ® J^)- 
// 77' is not of full rank, then 77' is a product operator if and 
only ifr\ is a product operator. 

Proof. Suppose that 77' is not full rank, and is a product opera- 
tor. As it is positive semidefinite, 77' must either have the form 
|a)(a| ® 77" or 77" ® \ot)(a\ for some \a) S Tf.2- In particular, 
there must exist a state vector I7) 6 H2 such that one of 

((7I ® 1)t?' (I7)® 1) =0 or (A3a) 
(1<8<7|K(1®|7» =0 (A3b) 



14 



holds. Decompose r\ in its spectral decomposition, 



then have 



V = ^2^k\<t>k){<j)k 

k 



forAfc > 0. Suppose that ((7I ® l)rf( 
|T) = U I7), we have 



(A4) 



= 0: if we let 



= ®1)(1®»7)(17|7)®1) 

= ((r|®i)^A*i®|^)(0 fc M(|r)®i) 

= ^AfeT^, (A5) 



where we define T k = ({T| ® ® |0fe)). By Lemma 1 
each operator Tfc is zero only if both |T) and \(f>k) are both 
product vectors; we may then decompose |T) = \a) |r) and 
\4>k) = \t') \4>'k)^ where we require (t|t') = for all k. We 



\r'){r'\ 



L k 



(A6) 



On the other hand, if (l ® (7] ® I7) ) = 0, we obtain 

= (t/ f ®<7|)r?(Z7® | 7 )) 
= ^A fc l7t(l®|^))(l®<^|)cr, 



(A7) 



for single-spin states \cf)' k ) = (l ® (7I ) l^). We then require 
t/t (l ® |</>' fc ) ) =0 for each k; because U cannot be decom- 
posed as U' ® \u) for any state \u) G 7^2, this implies that the 
vectors |<^) themselves are zero. Thus, \<j>k) — \ctk)\/3) for 
some states \ctk) E H2 , and where ("f\f3) = 0. We then have 



£ 



(A8) 



In either case, 77' is a product operator only if 77 is a product 
operator; the converse holds trivially. ■ 



